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HE: X 射线 暴 振荡 是 在 小 质量 X 射线 双星 的 I 型 X 射线 暴 中 探测 到 的 周期 性 现象 ， 振 荡 的 周 
期 与 中 子 星 的 自转 周期 接近 。 近 些 年 来 ，RXTE 和 NICER 等 高 能 探测 器 在 30 个 小 质量 X 射线 
双星 中 探测 到 了 爆发 振荡 现象 。X 射线 暴 振荡 对 研究 致密 物质 状态 、 中 子 星 表 面 强 引力 和 强 磁场 
行为 等 都 有 着 重要 的 意义 。 文 章 介 绍 了 I 型 X 射线 暴 和 爆发 振荡 的 搜索 方法 、 爆 发 振荡 的 观测 
Y 结果 、 爆 发 振荡 的 热 斑 模型 和 表面 模式 振荡 模型 ， 未 来 ，X 射线 暴 振荡 将 用 于 测量 中 子 星 质量 、 


半径 和 自转 。 
x 键 dg: X HE X 射线 暴 振荡 ， 小 质量 X 射线 双星 ; 脉冲 星 
中 图 分 类 号 : P142.6 文献 标识 码 : A 

1 引 言 


小 质量 义 射线 双星 (low-mass X-ray binary, LMXB) 由 致密 星 (中 子 星 或 黑洞 ) 和 一 颗 小 
质量 主 序 星 (M < Mo) 组成， 其 X 射线 辐射 来 自 于 致密 星 的 吸 积 过 程 。 如 果 致 密 星 为 中 子 
星 ， 吸 积 物质 累积 在 高 温 、 高 压 和 强 引 力 场 的 星体 表面 ， 可 触发 不 稳定 的 核 燃 烧 ， 观 测 上 表 
MIIR HRR, Strohmayer 等 人 "借助 罗 西 X 射线 时 变 探索 者 (Rossi X-ray Timing 
Explorer, RXTE) 高 时 间 分 辨 率 和 大 有 效 面积 的 特性 ， 首 次 探测 到 LMXB 4U 1728—34 I78 
X 射线 暴 中 频率 为 363 Hz 的 周期 性 信号 ， 即 X 射线 暴 振 荡 现 象 (X-ray burst oscillation, 
XBO)。2003 年 Chakrabarty 等 人 发 现 了 吸 积 毫秒 脉冲 星 SAX J1808.4—3658 在 I 型 义 射 
线 暴 中 的 周期 性 信号 ， 且 在 暴 尾 部 ，X 射线 暴 振荡 频率 与 自转 频率 只 偏差 约 6x10-? Hz; E 
年 ，Strohmayer 等 人 "在 另 一 颗 吸 积 毫秒 脉冲 星 XTE J1814—338 中 发 现 了 类 似 的 现象 。 这 
些 观测 证 实 了 X 射线 暴 振 荡 源 为 核 供 能 脉冲 星 ， 即 脉冲 辐射 能 量 来 自 星体 表面 核 燃 烧 。 
通过 RXTE, Swift 和 NICER 卫星 的 观测 数据 ， 发 现 LMXB 源 中 约 2096 89 L78 X 5512 
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展 荡 现象 ， 振 荡 的 频率 范围 约 为 11 ~ 620 Hz， 振 荡 的 均 方 根 振幅 在 


596 ~ 20% 之 间 ， 在 有 些 源 中 ， 均 方 根据 幅 可 接近 50%。 在 吸 积 毫 秒 脉冲 星 中 ， 观 测 到 的 


XBO 频率 接近 星体 的 自转 频率 ， 差 异 在 几 赫 兹 之 内 。 在 I 型 X 射线 暴 的 上 升 和 衰减 尾部 ， 


XBO 的 频率 有 几 赫 效 的 漂移 ， 通 常 在 


射线 暴 的 尾部 探测 到 。 

Strohmayer 等 人 ”认为 X 射线 暴 振 荡 的 信号 来 自 于 中 子 星 的 自转 ， 即 在 I 型 X 射线 暴 
爆发 时 ， 可 能 在 中 子 星 表面 形成 热 班 ， 导 致 不 均匀 的 温度 分 布 ， 由 星体 自转 在 X 射线 能 段 
产生 周期 性 信号 ， 但 I 型 X 射 线 暴 的 爆发 过 程 中 热 斑 如 何 形成 和 演化 等 问题 尚未 解决 。 表 


本 文 分 别 从 I 


ili 
同 高 度 形成 温度 梯度 ， 通 过 星体 自转 调和 
现象 ， 有 待 进一步 的 完善 。 


昆 部 渐 近 中 子 星 的 旋转 频率 ， 且 XBO 更 容易 在 I 型 又 


面 模式 模型 (如 Rossby 模式 ) 认为 1 型 X 射线 暴 可 激发 中 子 星海 洋 的 振荡 ， 在 星体 表面 不 
前 ， 形 成 XBO。 这 些 模型 只 能 解释 部 分 观测 的 XBO 


AI K WRT, X 射线 暴 振 荡 的 观测 以 及 振荡 理论 分 别 进行 介绍 。 


2 IAN X 射线 暴 


I7] X 射线 暴 


是 发 生 在 中 子 星 LMXB 中 的 不 稳定 核 燃 烧 ， 主 要 有 HANS He KR C 
燃烧 ， 其 中 C 燃烧 在 观测 上 为 超 暴 ， 由 Grindlay 等 人 ”在 4U 1820-30 中 首次 发 现 。 典 型 的 
I X 射 线 暴 的 光 变 曲线 旦 快速 上 升 和 指数 下 降 形状 ， 上 升 时 间 约 1 ~ 10 s ( 见 图 四 ， 持 
续 时 间 10 ~ 100 s， 爆 发 间隔 时 间 为 几 小 时 到 几 天 ， 能 谱 一 般 由 黑体 谱 拟 合 ， 特 征 温度 为 


1 ~ 3 keV， 辐 射 的 总 能 量 1032 ~ 1033 J。 目 前 在 115 个 中 子 星 LMXB 源 中 观测 到 了 I 型 X 


HREN, 


iE: a) 是 1996 4 


E 2 16 € 4U 1728-34 暴 的 光 变 曲 


Phttps://personal. 


sron.nl/rjeanz/bursterlist.html 


345 350 355 360 365 370 375 380 
f/Hz 
b) 
2X, b) 是 功率 密度 谱 中 的 XBO 信号 。 
1 I 型 义 射 线 暴 的 光 变 曲线 及 其 功率 密度 谱 " 


C— 
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TA X 射线 暴 的 核 燃 烧 过 程 敏 感 地 依赖 于 中 子 星 的 局 部 吸 积 率 ， 定 义 局 部 爱 丁 顿 吸 积 
率 为 ， 


2myc — . 15x10 g.em 2.s-1 (* m) (1) 
1 + X)Roqu 1-X R ' 

Ru. my 为 质子 质量 ，c 为 光速 ，X 为 大 气 中 的 耳 丰 度 ，om 为 汤姆 孙 散 射 截面 ， 尺 为 中 子 
星 半 径 。 对 于 太阳 丰 度 的 大 气 (X = 0.73) 且 半 径 为 10 km 的 中 子 星 ， 其 局 部 爱 丁 顿 吸 积 率 
为 8.8x104g:.cm-2.Ss-1。 

2.1  H 燃烧 

中 子 星 表面 温度 较 低 时 ，H 燃烧 由 冷 CNO 循环 主导 ， 少 数 情况 会 发 生 热 CNO 循环 。 
而 中 子 星 表面 温度 足够 较 高 时 ，H 燃烧 主要 是 rp 过 程 。 冷 CNO 循环 过 程 如 下 : 


12C+p 一 SN+TY 
13N — 130. Let Lp, 
SC+p > “N+Y 
MN+p 一 150 
150 > MN pet pu, 


BSN+p 一 ?Cca 


循环 过 程 中 ，C, N, O 元 素 捕获 质子 ， 并 在 反应 的 最 后 一 步 释 放 a 粒子 ，C, N, O 的 生成 速 
率 与 消耗 速率 相同 ， 是 循环 过 程 的 催化 剂 。 冷 CNO 的 反应 速率 依赖 于 捕获 质子 的 速率 ， 尤 
其 是 MN(p, VPO 反应 。 当 吸 积 率 低 于 0.017rgaa，H 经 冷 CNO 循环 过 程 燃烧 ， 反 应 速率 
非常 快 ， 这 个 过 程 释放 的 能 量 与 T1 ( 为 环境 温度 ) 成 正比 ， 反 应 过 程 不 稳定 ， 同 时 会 触发 
He 燃烧 ， 观 测 上 表现 为 H/He 混合 燃烧 工 型 X 射线 暴 。 
2.2 He 燃烧 

He 丰 度 较 高 的 壳 层 可 以 在 稳定 的 HH 燃烧 ( 吸 积 率 0.01rngaa ~ O.1rgaa) 的 壳 层 之 下 形 
成 ， 或 者 通过 吸 积 He 白矮星 、 大 气 成 分 为 He 的 主 序 星 形成 。 在 He 壳 层 中 ， 当 吸 积 率 为 
0.1mgaa ~ 1mrpaa， 环 境 温度 较 高 (T > 108 K)， 简 并 的 电子 密度 达到 105 g/cm? 时 ， 会 触 
发 不 稳定 的 He 燃烧 (反应 速率 oc T39， 此 过 程 极其 不 稳定 )， 此 时 He 燃烧 过 程 是 a 十 a 十 a 
f > C+y。 这 个 过 程 分 两 个 阶段 ， 第 一 阶段 是 两 个 a 粒子 结合 成 8Be，8Be 的 寿命 只 有 
© 2.6 x 10716 s， 此 时 会 源源 不 断 地 形成 8Be， 从 而 有 足够 的 5Be 粒子 与 a 结合 成 C 的 激发 态 
120+, 120* 可 以 衰变 成 8Be 并 释放 a 粒子 ， 或 者 辐射 y 光子 变 成 C 基态 ， 这 是 3o 过 程 的 
第 二 个 阶段 。 观 测 上 表现 为 He 燃烧 I 型 X 射线 暴 ， 与 H/He 混合 燃烧 I 型 X 射线 暴 相 比 ， 
其 光 变 曲线 的 上 升 下 降 时 间 更 短 。 
2.3 C 燃烧 

在 很 少 的 几 个 源 中 (例如 4U 1820-30)， 观 测 到 一 类 特殊 的 X 射线 暴 ， 也 称 为 超 暴 
(superburst)， 与 上 述 两 类 I 型 X 射线 暴 相 比 ， 超 暴 释放 的 总 能 量 高 约 109 倍 ， 持 续 时 间 长 
约 108 倍 ， 发 生 率 约 为 -5， 这 类 暴 可 能 是 中 子 星 大 气 中 更 底层 的 C 壳 层 的 燃烧 。 目 前 ， 
超 暴 的 发 生 条 件 尚 不 清楚 ， 可 能 需要 额外 的 加 热机 制 (例如 ， 中 子 星 壳 层 的 超 密 态 核反应 )， 


meqd = ( 
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以 克服 C REMEC R. 


3 X 射线 其 振荡 的 


对 XBO 的 探测 ， 


死 时 间 效 应 不 显著 等 特性 。 


Z E E bb 
需要 高 能 


E 1996 年 RXTE 首次 发 


观测 


又 射线 望远镜 


30 个 源 中 发 现 了 这 类 现象 ，XBO 频率 见 表 1， 其 中 确定 暴 振荡 频率 的 源 按 脉 冲 星 和 非 脉冲 


星 分 类 ，XBO 频率 有 待 进 
MI NICER 卫星 的 观测 中 搜索 到 了 XBO fs. FH, > 


号 的 概率 估计 和 XX 射线 望 


远 镜 对 XBO 的 观测 。 


x 


步 观测 验证 的 源 单独 分 出 。 


39 3$ 


了 高 时 间 分 辨 率 、 大 有 效 面 积 、 光 子 堆积 和 
现 中 子 星 LMXB 的 XBO 起 ， 目 前 共 在 


除 RXTE 外 ， 人 们 还 从 Swift/BAT 
各 分 别 介绍 XBO 信号 的 搜索 方法 、 信 


表 1 ”探测 到 XBO ME" 7T Hz 
爆发 源 振荡 频率 ”自转 频率 探测 器 参考 文献 
脉冲 星 

SAX J1808.4—3658 401 401 RXTE, NICER 四 , c3 
IGR J17498—2921 401 401 RXTE E3, a] 
XTE J1814—338 314.4 314 RXTE 回 
IGR J17511—3057 245 245 RXTE [ES 
Aql X-1 549 550.27 RXTE Id 
HETE J1900.1—2455 376.25 377 RXTE, Swift I 
IGR J18245—2452 254.4 254 Swift [E] 
非 脉 冲 星 
4U 1608 一 522 620 RXTE Il 
SAX J1750.8—2900 601 RXTE [g, uj 
GRS 1741.9—2853 589 RXTE IZ T] 
4U 1636—536 581 RXTE [22]. 世相 (在 一 个 超 暴 里 ) 
X 1658—298 567 RXTE [af 
KS 1731—260 524 RXTE 23, Zu] 
4U 1728—34 363 RXTE, NICER m, zz 
4U 1702—429 329 RXTE [Z3] 
IGR J17191—2821 294 RXTE [29] 
SAX J1810.8—2609 532 RXTE [Ed 
EXO 0748—676 552a RXTE BI, g3 
IGR J17480—2446 11 RXTE [33] 
XBO 频 率 有 待 证 实 
XTE J1739—285 11225 RXTE Ba 
GS 1826—24 611° RXTE Ba 
A 1744—361 529% RXTE Eo 
SAX J1748.9—2021 409.7° 442.36 RXTE Bz 
4U 0614+09 415f RXTE, Swift B3] 
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( 续 表 ) 
爆发 源 振荡 频率 ”自转 频率 探测 器 参考 文献 
MXB 1730-335 306.58 RXTE [39] 
4U 1916—053 270^ RXTE zol 
XB 1254 一 690 95! RXTE g 
ik: 上 标 “a” 表 示 2004 ££ Villarreal 和 Strohmayer 测 得 的 45 Hz 的 振荡 信号 ， 而 2010 年 Galloway 


等 人 四 又 测 得 552 Hz 的 振荡 信号 ， 并 且 Galloway 用 同样 搜索 方法 去 搜索 2004 年 的 129 个 爆发 样本 ， 并 
没有 发 现 45 Hz 的 振荡 频率 。552 Hz 的 振荡 特性 与 其 他 源 中 所 观测 到 的 暴 振荡 频率 、 振 幅 以 及 动态 功率 谱 相 
似 ， 他 们 认为 552 Hz 更 接近 EXO 0748-676 中 子 星 的 自 旋 频 率 。 

上 标 “b” 表 示 Kaaret 等 人 四 在 6 个 义 射线 暴 中 探测 到 1122 Hz 的 振荡 信号 。 但 在 同样 的 数据 中 ， 其 
他 人 并 没有 找到 这 个 振荡 信号 ，Kaaret 等 人 得 到 振荡 信号 可 能 是 噪声 模型 不 够 准确 。 

上 标 “c” 表 示 Thompson 等 人 四 在 3 个 暴 的 蕉 加 谱 中 发 现 611 Hz 的 振荡 信号 。Bilous 和 Watts” 
在 77 个 暴 中 搜索 时 ， 没 有 在 611 Hz 附近 搜索 到 振荡 信号 。 
上 标 “d” 表 示 Bhattacharyya 等 人 E RIL 529 Hz 的 振荡 信号 ， 但 在 后 面 对 真 实 信 号 概率 进行 估计 时 
RI, RA 196 的 模拟 具有 振荡 信号 。 

上 标 “e” 表 示 Kaaret 等 人 四 发 现 409.7 Hz 的 振荡 信号 ， 但 是 Altamirano 等 人 四 在 442.36 Hz 处 
发 现 了 强 脉冲 信号 ， 所 以 认为 Kaaret 等 人 发 现 的 振荡 信号 不 够 准确 。 

上 标 “f” 表 示 Strohmayer 等 人 四 在 暴 尾部 发 现 415 Hz 的 振荡 信号 ， 而 RXTE 只 从 一 个 暴 中 的 一 个 
时 间 窗 口中 搜索 到 爆发 振荡 。 
上 标 “g” 表 示 Fox 等 人 四 在 31 个 暴 关 加 的 功率 谱 上 发 现 306.5 Hz 的 爆发 振荡 信号 。 但 是 Bilous 和 
Watts"! 在 57 个 爆发 数据 中 没有 在 306.5 Hz 附近 搜索 到 任何 振荡 信号 。 
上 标 “h” 表 示 Galloway 等 人 四 从 6 个 暴 中 搜索 到 270 Hz 的 暴 振荡 信号 ， 但 是 没有 考虑 到 时 频 窗 口 数 
量 的 重要 性 ， 并 且 归 一 化 之 后 Pu 的 值 降低 很 多 四 。 
上 标 “i” 表 示 Bhattacharyya 搜索 到 95 Hz 的 暴 振荡 信号 ， 但 是 Bilous 和 Watts” 认为 不 排 
除 是 噪声 产生 的 结果 。 


3. R XBO 的 方法 

搜索 周期 性 的 XBO 信号 ， 可 通过 快速 傅 里 时 变换 (fast fourier transformation, FFT) 
或 22 统计 方法 ”， 它 们 分 别针 对 并 道 数据 及 事例 数据 加 . 

对 于 FFT， 将 总 持续 时 间 工 的 光 变 曲线 均匀 划分 为 NN 个 并 道 ， 且 N = 2" (m 是 整数 )， 
时 间 序 列 记 为 zs (k = 0,1,2,3,- ,入 一 1)， 表 示 第 个 并 道中 的 光子 数 ，Leahy 归 一 化 功 
率 密度 谱 P; 为 ” : 


N N 


2 2 
2 
po AUS (> Tk COS 22 十 (> Tk Sin znt : (2) 


k=1 k=1 


Ep, y 是 传 里 叶 频 率 (v; = j/T, j = 0,1,2,3,...,N/2) Mo 为 总 光子 数 。 功 率 密度 谱 
的 最 小 频率 为 Vmin = 1/T， 最 大 频率 ( 即 奈 斯 奎 特 频 率 ) 为 vma = N/(2T)， 频 率 间隔 为 
Av =1/T. 


© 


jn 


友人 例 数据 记录 了 每 个 光子 的 到 达 时 间 ， 且 可 转换 为 并 道 数据 ;并 道 数据 记录 每 一 个 时 间 间 隔 中 的 光子 计数 率 。 
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Z2 方法 与 FFT 计算 的 功率 密度 谱 的 过 程 相似 ， 它 的 定义 为 ”: 


Ny Ny 2 
E = >， lp cos tos) 一 (X sin 2 | , (3) 
tot k=1 | \j=1 j=1 


其 中 ，n 为 谐 波 数 ，9; 为 光子 相位 : 


2 


i cn f ' v(t)dt. , (4) 


其 中 ，to 是 参考 时 间 ，t; 是 光子 相对 于 to 的 到 达 时 间 ，v(t) 是 频率 模型 ， 当 v(t) 不 随时 间 
变化 时 ，9; = 27vt;j。 

除了 振荡 频率 ， 功 率 谱 还 包含 有 关 脉 冲 幅度 的 信息 。 对 于 Leahy 归 一 化 的 功率 谱 ， 均 
方 根 分 数 (root mean square, rms) 振幅 定义 为 : 


P. N. 
A 5 ii 
id N, (x = z) ? (5) 


括号 项 用 来 校正 背景 因素 影响 ，y 是 总 光子 数 ，B 是 背景 光子 数 四 ， 已 , 是 信号 功率 ， 如 果 
使 用 Z2 方法 ， 则 将 公式 (B) 中 的 P. 蔡 换 为 Zio 

与 FFT 相 比 ，2% 方法 能 更 精准 地 测 得 振荡 频率 ， 但 计算 时 间 更 长 。 通 常 ， 对 于 I 型 X 
射线 暴 的 事例 数据 ， 先 转换 为 并 道 格 式 ， 然 后 使 用 FEFT 方法 搜索 数据 中 是 否 存 在 XBO 信 
号 ， 如 果 存 在 (或 有 疑似 信号 )， 可 进一步 使 用 2 方法 更 准确 地 测量 振荡 频率 。 
3.2 ”如 何 估计 真实 信号 概率 

用 FFT 或 者 Z2 统计 方法 找到 候选 的 XBO 信号 之 后 ， 需 要 判断 信号 真实 的 概率 ， 男 
外 ， 也 需要 由 测量 的 XBO 信号 给 出 真实 信号 的 概率 分 布 。 
3.2.1 估计 快速 傅 里 叶 变 换 的 真实 信号 概率 

由 了 FT 方法 计算 I 型 X 射线 暴 的 功率 密度 谱 ， 低 频 包含 红 噪声 和 泊 松 噪声 ， 而 高 频 部 
分 主要 由 泊 松 噪声 组 成 。 泊 松 噪声 产生 的 功率 密度 具有 两 个 自由 度 的 x? 分 布 ”: 

2(R2 + I 
AT ENS ， M 

Jt, m, L, 分 别 为 噪声 傅 里 叶 变 换 系 数 的 实 部 和 虚 部 ， 根 据 Groth” 解析 表达 式 ， 测 量 
功率 Py, 和 信号 功率 P, 的 概率 分 布 : 


(k-1)/2 
1 LB Fr 
px(Pn, Ps) = gu ( 2 j (3 ) X Iu (EF. ) (7) 


其 中 ,了 工 是 第 一 类 的 修正 贝 塞 尔 函 数 ， 而 天 是 谐 波 数 (一 般 取 1)。 对 于 仅 有 噪声 的 情况 ， 
PP = 0, 有 == 1， 公 式 (@) 为 具有 两 个 自由 度 的 X? 分布 的 概率 密度 函数 (probability density 


“一般 取 暴 前 的 光子 平均 计数 率 乘 以 计算 功率 密度 谱 所 对 应 的 时 间 段 宽度 ， 作 为 背景 光子 数 
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function, pdf)， 等 于 
1 =P a 
pi(Pn,0) = 3 P (= ) ; (8) 


XF k > 1, pr(Pa,0) 满足 自由 度 为 2k 的 x? 分 布 。 
3.22 ”估计 Z2 统计 的 真实 信号 概率 

Z2 统计 判断 由 噪声 产生 信号 几率 小 于 1%， 则 视 为 真实 信号 ， 以 Ootes 等 人 在 2017 
年 所 处 理 的 RXTE 爆发 数据 方法 为 例 ， 把 检测 范围 定 为 测 得 大 概 的 振荡 频率 前 后 5 Hz. "ul 
回 所 示 ， 将 试验 的 次 数 定 义 为 寻找 信号 的 时 间 和 频率 (时 频 ) 窗口 的 总 数 ，N = Ni x Ny, 
Ny, 为 时 间 窗 口 的 数量 ，NN, 为 频率 窗口 的 数量 。 


Vi 十 5 


时 间 窗 口 


0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 
T/s 


ik: 搜索 暴 振荡 频率 而 创建 的 时 频 窗口 可 视 化 ， 爆 发 的 前 16 s 为 时 间 窗 口 ， 每 个 时 间 窗 口 对 应 10 个 频率 窗口 是 
已 知 频率 的 前 后 5 Hz， 每 个 时 间 窗 口 都 有 5000 个 计数 ， 所 以 时 间 窗 口 在 时 间 坐 标 上 的 分 布 不 均匀 ， 颜 色 刀 
不 同 标准 。 蓝 色 对 应 于 标准 1 (单个 时 间 窗 口 进 行 检测 )， 红 色 对 应 于 标准 2 (对 第 一 个 时 间 窗 口 进行 检测 )， 绿 
对 应 于 标准 3 ( 取 连 续 两 个 时 频 窗 


Em 


KI 


图 2 ”时 频 窗口 可 视 化 下 


假设 噪声 为 泊 松 分 布 ， 由 噪声 产生 噪声 分 布 0 的 测量 信号 概率 Poiso 73: 


Poise = NO6(1 — 6AT , (9) 


根据 Poise < 1%，Muno 5 A" 所 定义 的 三 个 标准 来 判断 重要 信号 : 

(1) 将 N= N, x N, = 16 x 10 = 160 KEI Paoise < 1% RAR ( 团 ， 可 以 得 到 噪声 
分 布 5 LT x 1075. 

(2) 计算 暴 上 升 阶段 第 1 s， 即 Ni = 1, Paise < 1% 时 由 公式 (B) 得 到 噪声 分 布 
6 < 1073, 

(3) 一 个 时 频 窗口 周围 有 两 个 相 邻 的 频率 窗口 以 及 一 个 时 间 窗 口 ， 计 算 所 选 时 频 窗 口 和 
其 中 一 个 相 邻 窗口 的 噪声 概率 ; 两 个 窗口 单独 包含 噪声 的 概率 由 两 个 窗口 独立 的 噪声 概率 
的 乘积 给 出 ，Probls = Pi(N1,01) x Po(N2, 02); AETH 3 个 相 邻 窗口 计算 噪声 概率 ， 所 以 
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其 中 一 个 和 N, = 3， 可 近似 表示 为 : 
Prob; > & 3N? N, 040; ; (10) 
M Proba 2 <10? 时， 将 其 代入 式 (四 )， 则 有 6162 < 1.3 x 1079, 
这 些 标准 以 及 Z2 统计 满足 自由 度 为 2n 的 x? 分 布 (n XS BC 1), AA Z2 统计 中 
的 功率 是 否 为 真实 信号 ， 对 于 给 定 的 信号 功率 Zo MEAS Zm 的 条 件 概率 分 布 py WEA 
XU. 
1 Ww x gm y cue 
Dk (Zm|Zs) = 2 exp | ( m s) (2) X Iki (v ZZ.) ; (11) 
其 中 , k EWEA (k= 1)， 而 Doa 是 第 一 类 修正 的 Bessel 函数 。 给 定 2Z,， 通 过 相关 


的 累积 分 布 函数 获得 介 于 0 和 Zm 之 间 的 测量 功率 的 概率 是 : 


Zm Zs) oo kn— 2 
fr (ZmlZs) = exp | H | x p PE ME (12) 
在 给 定 Zm 的 情况 下 ， 真 实 信号 功率 在 0 到 Z, 间 的 概率 是 : 
fk (Zs|Zm) = 1— fr (Zm|Zs) (13) 
如 果 在 整个 爆发 中 均 未 找到 振荡 信号 ， 则 选择 最 大 的 功率 作为 信号 的 上 限 。 
3.238 蒙特 卡 洛 模拟 
蒙特 卡 洛 模拟 方法 可 用 来 评估 无 信号 时 噪声 功率 的 分 布 ， 且 能 精确 模拟 实际 观测 的 数 
据 ， 如 设 定 的 时 间 窗 口 等 。 首 先 选 择 合 适 的 模型 拟 合 I 型 X 射线 暴 的 光 变 曲线 ， 然 后 以 此 


为 基础 生成 大 量具 有 泊 松 计数 统计 信息 但 没有 任何 周期 性 信号 的 光 变 曲线 样本 。 从 这 些 模 
拟 出 来 的 光 曲 线 中 获取 功率 谱 ， 即 可 得 出 仅 由 噪声 构成 的 功率 密度 谱 。 在 测 得 的 振荡 信和 号 附 


近 统 计 功 率 分布 ， 此 分 布 即 为 噪声 
3.8 XBO 的 动态 功率 密度 谱 
Strohmayer 等 人 ”将 整个 I 型 X 射线 暴 的 光 变 曲线 应 用 FFT 方法 ， 发 现 功率 密度 谱 中 
结构 ( 见 图 加， 这 是 因为 在 整个 爆发 期 间 XBO 的 振荡 频率 随时 间 演 化 。 为 了 有 
效 地 搜索 和 追踪 XBO 信和 号， 一般 用 动态 功率 密度 谱 方法 分 析 工 型 X 射线 暴 。 
图 加 给 出 了 产生 动态 功率 谱 的 示意 图 ， 图 Bla) 为 1 型 X 射线 暴 的 光 变 曲线 ， 可 将 光 变 曲 
线 划 分 为 多 段 ， 这 里 以 三 段 为 例 ，Ato, Ati, At: 是 移动 窗口 的 时 间 长 度 ， 相 邻 的 窗口 之 间 
FEZ, KE TD, Dn Mn 为 窗口 的 中 间 时 刻 。 图 苞 ) 分 别 计算 每 段 光 变 曲线 的 功率 密度 谱 ， 
蓝 色 代表 低 功 率 ， 紫 色 代 表 中 功率 ， 红 色 代 表 高 功率 。 图 多) ES Bo) 的 功率 密度 谱 和 相对 
应 的 中 间 时 刻 结合 ， 给 出 时 间 - 功 率 -频率 三 维 图 。 最 后 在 图 名) E ge) 功率 的 等 高 线 投 
影 在 光 变 曲线 上 ， 得 到 动态 功率 密度 谱 ， 可 简洁 地 展示 XBO 出 现 的 时 间 、XBO 信和 号 的 频 
率 和 强度 随时 间 的 演化 等 信息 。 


分 布 ， 可 判断 出 测量 到 的 XBO 信号 来 自 于 噪声 的 概率 。 


出 现 了 多 峰 
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li) 


: 
E 


时 间 


Hs 图 a), b), c), d) 表示 动态 功率 谱 的 搜索 过 程 。 首 先 选 取 移动 时 间 窗口 ， 计 算 每 段 光 变 曲线 的 功率 密度 谱 ， 
之 后 将 光 变 曲线 和 功率 密度 谱 结合 到 一 起 ， 最 终 得 到 动态 功率 谱 。 


图 3 动态 功率 密度 谱 


3.4 RXTE 对 XBO 的 观测 

RXTE 卫星 于 1995 年 12 月 30 日 发 射 升 空 ， 在 长 达 16 年 的 运行 后 ， 于 2012 年 1 月 5 
日 正式 退役 ， 该 卫星 是 到 目前 为 止 研究 XBO 最 重要 和 最 成 功 的 高 能 时 变 探测 器 。RXTE E 
星 采 用 事例 格式 记录 高 时 间 分 辩 率 数据 ，2019 4E, Bilous 和 Watts" 用 FFT 方法 系统 搜索 
和 分 析 了 所 有 RXTE 观测 的 I 型 X 射线 暴 的 XBO 信号 ， 搜 索 了 所 有 57 个 中 子 星 LMXB 
源 的 2118 次 爆发 ， 列 出 了 爆发 峰值 的 信 噪 比 ， 提 取 了 2 ~ 60 keV 能 量 范围 内 的 光 变 曲线 ， 
搜索 时 间 窗 口 为 2 s， 移 动 步 长 为 0.5 s， 提 供 了 暴 振荡 频率 以 及 在 爆发 中 出 现 的 位 置 等 信 
息 ， 并 计算 了 振荡 的 分 数 振幅 ， 对 已 经 确定 振荡 频率 的 源 进行 搜索 ， 确 定 了 已 知 的 17 个 
XBO 频率 与 之 前 测 的 频率 或 者 脉冲 星 自传 频率 相等 ， 对 之 前 一 些 论文 中 所 搜索 到 的 爆发 振 
荡 但 还 有 待 检验 的 源 进 行 检 验 ， 其 中 一 些 源 并 没有 搜索 到 暴 振荡 信号 ， 另 一 些 源 搜索 到 的 
振荡 信号 过 于 微弱 ， 不 排除 是 噪声 信号 (关于 XBO 频率 有 待 证 实 的 源 的 具体 信息 ， 可 见 表 
15 还 发 现 了 SAX J1810.8—2609 在 531 Hz 的 爆发 振荡 信号 以 及 IGR J17473—2721 频率 在 
600 Hz 左右 的 振荡 信号 ， 该 源 不 同 暴 的 振荡 频率 之 间 相 差 在 3 Hz 以 内 ， 很 可 能 是 一 个 新 的 
XBO 源 ， 其 准确 振荡 频率 还 需要 进一步 检验 。 
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而 Galloway 等 人 主要 用 Z2 方法 搜索 ”， 其 中 包括 85 个 中 子 星 LMXB 源 ，7083 次 爆 
发 ， 主 要 取 3 ~ 25 keV 能 量 范围 内 的 光 变 曲线 进行 搜索 ， 最 后 在 16 个 不 同 的 源 中 搜索 到 
244 次 爆发 振荡 ， 提 供 了 爆发 振荡 信号 真实 性 的 概率 以 及 频率 随时 间 演 化 、 检 测 到 最 强 信和 号 


的 爆发 相位 等 信息 。 


图 团 是 RXTE 观测 到 的 典型 的 XBO 动态 功率 密度 谱 ， 其 中 图 Ba) SAX J1808.4—3658 


MA Ab) XTE J1814—338 是 两 个 持续 吸 积 脉冲 星 ， 图 欧 ) HETE J1900.1—2455 VX BIM) 


Aql X—-1 为 间 欺 吸 积 脉 冲 星 ， 而 图 区 ) 4U 


404 


T 


SAX J1808.4-3658 


f/Hz 


T/s 
e) 


1636—536. K| mf) 4U 1728—34 都 是 非 脉 冲 星 。 


4 


V XTE J1814-338 i 


m 
F | 


计数 率 /103.s-1 


ik: RXTE 对 典型 的 LMXB 的 XBO 动态 功率 密度 谱 的 观测 。 绿 色 线 为 光 变 曲线 ， 蓝 色 为 功率 谱 的 等 高 线 ， 


虚线 为 脉冲 星 的 自转 频率 。 其 中 a) 和 b) 是 两 个 持续 的 吸 积 毫秒 脉冲 星 ，c) 和 d) 分 别 为 间歇 吸 积 脉 冲 星 ; e) 和 


£) 是 为 两 个 非 脉冲 星 。 


图 4 ”不 同类 型 源 的 动态 功率 密度 谱 ” 


研究 发 现 SAX J1808.4—3658 在 暴 开 始 时 就 出 现 振荡 信号 ， 并 且 随 着 暴 进行 向 上 或 者 向 
下 有 几 Hz 的 漂移 ， 在 暴 峰 值 振荡 突然 消失 ， 之 后 振荡 会 在 稍 高 或 稍 低频 率 处 重新 出 现 并 保 
fa: XTE J1814—338 的 振荡 频率 变化 在 一 个 傅 里 叶 频 率 分 辨 率 之 内 ， 振 荡 信 号 也 不 会 


在 暴 峰 值 期 间 消 失 ， 而 且 经 常 出 现在 暴 的 
相近 ， 在 几 Hz 之 内 。HETE J1900.1—245 


其 余 频 率 比较 低 的 振荡 信号 会 持续 一 段 时 间 ; Aql X-1 和 HETE J1900.1—2455 等 间歇 吸 积 


尾部 ， 持 续 吸 积 脉 冲 星 振荡 频率 与 自转 频率 非常 
5 只 有 一 个 时 间 窗 口 检测 到 爆发 振荡 信号 的 峰值 ， 


脉冲 星 的 爆发 振荡 频率 ， 在 整个 爆发 过 程 


FPF 有 缓慢 的 漂移 ， 渐 近 最 大 频率 且 比 自转 频率 低 约 
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1 Hz。4U 1636—536 爆发 振荡 信号 存在 于 暴 上 升 阶 段 ，4U 1728—34 的 振荡 信号 存在 于 整个 
爆发 过 程 中 ; 4U 1636 一 536 和 4U 1728—34 等 非 脉冲 星 所 观察 到 的 爆发 振荡 ， 随 着 爆发 的 进 
行 ， 与 间 葡 吸 积 脉 冲 星相 似 ， 频 率 约 1 Hz 向 上 漂移 ， 在 暴 的 尾部 渐 近 最 大 频率 。 

RXTE 还 从 中 子 星 LMXB 的 超 暴 中 发 现 了 XBO 现象 。Strohmayer 和 Markwardt"" 在 
2001 年 观测 到 的 4U 1636-53 的 超 暴 光 变 曲线 和 动态 功率 谱 见 图 回 在 动态 功率 谱 中 还 可 以 
看 到 明显 的 频率 漂移 ， 来 自 于 中 子 星 在 其 双星 轨道 上 的 多 普 勒 运动 。 超 暴 振荡 频率 比 一 般 I 
W X 射线 暴 的 振荡 频率 高 出 约 0.4 Hz， 比 正常 爆发 振荡 更 接近 中 子 星 的 真实 自转 频率 。 另 
外 ， 在 LMXB 4U 0614 上 +09 写 中 也 发 现 了 超 暴 的 爆发 振荡 信和 号 。 
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Hk: a) 是 2001 年 2 月 22 日 观测 到 的 4U 1636—53 的 超 暴 的 光 变 曲线 ， 时 间 分 辨 率 为 1 s。 垂 直 虚 线 表 示 在 检 
测 脉冲 的 800 s 间隔 观测 期 间 ， 源 被 地 球 遮 挡 ， 部 分 时 间 段 的 计数 率 为 0。b) 是 4U 1636—53 在 800 s 的 时 间 
间隔 内 的 动态 功率 密度 谱 ， 每 个 时 间 窗 口 为 64 s， 移 动 步 长 为 16 s。 


5 ”4U 1636-53 超 暴 的 光 变 曲 线 和 动态 功率 密度 谱 


© 35 “中子星 内 部 成 分 探索 器 (NICER) 探测 XBO 
I 2017 ££ 6 H 3 H, NASA 的 中 子 星 内 部 成 分 探索 器 (NICER) 发 射 成 功 ， 它 的 目标 是 通 
过 精确 测量 中 子 星 的 质量 和 半径 来 研究 中 子 星 的 物质 状态 和 内 部 结构 。NICER 提供 了 56 个 
准 直 X 射线 聚 光 器 阵列 ， 每 个 都 配 有 一 个 硅 漂移 探测 器 ， 这 些 探测 器 工作 在 0.2 ~ 12 keV 
能 段 ， 在 1 keV 处 的 能 量 分 辩 率 约 为 85 eV, TE 1.5 keV 处 的 有 效 面积 超过 2000 cm2， 时 间 
分 辨 率 小 于 300 ns， 显 著 地 优 于 RXTE. NICER 成 为 继 RXTE 之 后 研究 XBO 最 主要 和 最 
优秀 的 设备 。 
3.5.1 SAX J1808.4—3658 
Bult 等 人 在 2019 年 用 NICER 对 SAX J1808.4—3658 源 进行 观测 ， 发 现 了 预期 的 
401 Hz 频率 处 的 XBO， 且 没有 明显 的 漂移 。 他 们 还 给 出 了 XBO 的 脉冲 轮廓 和 持续 辐射 的 
脉冲 轮廓 ， 两 种 波形 如 图 日 所 示 ， 可 以 看 出 XBO 的 轮廓 和 幅度 与 持续 吸 积 脉冲 相似 ， 这 表 
明 两 者 都 源 于 中 子 星 表面 的 一 片 受 限制 区 域 的 燃烧 ( 热 班 ， 详 见 第 5E0 节 )， 但 是 XBO 领先 
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脉冲 34° 土 7? ，XBO 的 波形 与 能 量 的 相关 性 与 持续 吸 积 脉 冲 星 不 同 ， 一 般 在 低能 谱 测 得 暴 
发 振荡 信号 。 


波形 ， 


3.5.2 4U 
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ik: SAX J1808.4—3658 中 观察 到 的 暴 振荡 波形 与 其 持续 辐射 时 的 脉冲 波形 相 比较 ， 其 中 绿色 阴影 部 分 为 脉冲 


黑色 实 线 为 暴 振荡 波形 。 


Elo 暴 振荡 波形 与 脉冲 波形 相 比 较 图 


1728 一 34 


从 开始 工作 到 2019 年 2 HE, NICER 探测 到 了 4U 1728—34 的 7 次 I 型 X 射线 暴 ， 


Mahmoodifar EAT 发 现 其 中 的 3 次 爆发 有 XBO 现象 ， 且 都 出 现在 其 暴 衰 减 尾部 (如 图 口 


所 示 )， 频 率 非常 接近 ， 与 之 前 的 观测 结果 一 致 。 


比较 这 三 个 在 暴 尾部 的 振荡 ， 表 明 暴 在 高 能 带 中 可 以 更 强 ， 并 且 可 以 具有 不 同 的 幅度 和 
特性 。 其 中 的 两 个 振荡 只 在 高 能 段 (E > 6 keV) 中 被 检测 到 ， 并 且 它 们 具有 非常 大 的 分 数 均 


方 根 振幅 


ar 


(48: 9)96 和 (46 土 9)%。 热 斑 模 型 (第 5 节 ) 和 表面 模式 振荡 模型 (BEAT) 都 难 


以 解释 如 何在 爆发 衰减 尾部 出 现 的 大 振幅 的 振荡 。 如 果 中 子 星 表面 的 所 有 位 置 都 以 相同 的 
方式 加 热 和 冷却 ， 则 无 法 在 中 子 星 表面 产生 足够 高 的 温度 梯度 ; 而 冷却 尾 流 (cooling wake) 
模型 ” ( 详 见 E3 节 ) 表明 ， 如 需要 在 中 子 星 表 面 形成 很 高 的 温度 梯度 ， 则 要 中 子 星 表面 先 燃 
烧 部 分 比 后 燃烧 部 分 更 早 冷却 ， 在 燃烧 到 赤道 点 时 火焰 会 迅速 蔓延 至 整个 亦 道 ， 进 而 使 整个 
中 子 星 表面 的 温度 梯度 增加 ， 所 以 一 般 温 度 梯度 存在 于 赤道 以 及 附近 的 区 域 ”。 
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f/Hz 
计数 /s 


T/s 


7 NICER 观测 到 4U 1728—34 的 光 变 曲线 上 的 动态 功率 密度 谱 四 


4 WIRED X 射线 脉冲 星 与 XBO 的 关系 


吸 积 毫秒 脉冲 星 (accreting millisecond X-ray pulsars, AMXPs) 是 吸 积 供 能 (通过 吸 积 
伴星 的 物质 ， 将 引力 势能 转换 为 X 射线 辐射 ) 的 快速 转动 中 子 星 ， 正 是 在 AXMPs 中 观测 到 
了 XBO， 才 证 实 了 XBO 信号 来 自 于 中 子 星 自转 。 研 究 AXMPs 的 XBO 可 以 更 好 地 了 解 爆 
发 振荡 的 物理 机 制 。 这 类 源 的 自转 频率 可 以 被 独立 测量 ， 而 且 ， 它 们 是 可 以 用 来 研究 星体 磁 
场 影响 以 及 燃料 分 布 不 均匀 作用 的 唯一 一 类 源 。 

AMXPs 的 XBO 频率 与 其 自转 频率 非常 接近 ， 通 常 都 是 在 儿 Hz 以 内 ， 但 是 不 同 AMXP 
的 XBO 频率 与 自转 频率 的 间隔 存在 差异 。 例 如 XTE J1814—338 的 XBO 频率 非常 接近 自转 
频率 ; 而 SAX J1808.4—3658 和 IGR J17511—3057 的 暴 振 荡 在 暴 上 升 阶段 快速 增加 ， 之 后 
略微 超过 自转 频率 ， 最 后 在 暴 尾 部 会 稳定 在 自转 频率 为 10 mHz 左右 ”一 。AXMPs 与 非 肪 
冲 星 的 XBO 现象 也 有 差异 。 对 于 非 脉 冲 星 的 I 型 X 射线 暴 ， 当 源 处 于 软 态 (高 吸 积 率 ) 时 ， 
才 会 出 现 XBO 现象 ”"。 相 比 之 下 ， 即 使 在 硬 态 ( 低 吸 积 率 )，AXMPs 中 也 能 检测 到 XBO 
现象 AXMPSs 和 非 脉冲 星 的 暴 振荡 持续 时 间 也 存在 着 不 同 ， 对 于 SAX J1808.4 一 3658 和 
XTE 1814 一 338， 暴 振荡 现象 会 在 整个 暴 中 持续 ， 对 于 非 脉冲 星 ， 虽 然 有 时 会 在 整个 爆发 中 
检测 到 爆发 振荡 ， 但 是 大 多 暴 振荡 都 出 现在 暴 的 尾部 ” ， 这 点 和 间歇 吸 积 供 能 脉冲 星 一 样 。 
IGR J17511—3057 是 上 面 两 种 现象 的 过 渡 : 在 较 低 吸 积 率 下 ， 整 个 暴 中 都 可 以 检测 到 振荡 ， 
但 是 随 着 吸 积 率 的 增加 ， 爆 发 的 振荡 信号 会 从 暴 的 上 升 阶段 中 消失 ”。AXMPs 的 XBO 信 
号 比 间 葡 脉冲 星 和 非 脉冲 星 的 谐 波 含量 更 高 "”””， 并 且 它 们 振幅 与 能 量 的 关系 也 不 一 样 ， 
AXMP 的 XBO 振幅 随 能 量 的 上 升 而 下 降 ， 而 间歇 脉冲 星 和 非 脉冲 星 的 XBO 振幅 随 能 量 的 
E3Hfi EJF"^ 7, KBO 振幅 在 不 同 的 源 之 间 也 有 很 大 的 区 别 ， 大 多 数 情况 下 都 在 爆发 时 吸 
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职 供 能 脉冲 幅度 的 1% 到 10%; 但 是 对 于 过 渡 性 AMXP 射电 脉冲 星 IGR J18245—2452 
是 一 个 例外 : 对 于 检测 到 的 其 中 一 个 暴 ， 爆 发 的 振荡 幅度 远 高 于 爆发 由 吸 积 动 力 产 生 的 脉冲 
幅度 ”。 目 前 ， 对 造成 这 些 脉 冲 星 和 非 脉 冲 星 暴 振荡 之 间 特 性 差异 的 原因 尚 不 清楚 ， 动 态 
磁场 可 能 起 主要 作用 。 


5 XBO 理论 解释 


在 XBO 的 观测 中 发 现 ， 振 荡 信号 在 大 多 数 暴 中 频率 会 向 上 漂移 (频率 增加 )， 一 般 出 现 
在 暴 的 尾部 ， 在 一 些 暴 中 振荡 信号 有 着 很 高 的 振幅 ， 吸 积 脉冲 星 的 振荡 频率 与 脉冲 星 自转 
频率 的 偏差 在 几 Hz 之 内 ， 研 究 爆发 过 程 中 的 点 火 和 火焰 蔓延 方式 对 了 解 暴 振荡 的 物理 机 
制 起 着 重要 的 作用 。 目 前 提出 了 热 斑 模型 、 表 面 模式 振荡 模型 (如 图 加 和 冷却 尾 流 (cooling 
wake) 模型 来 解释 这 些 XBO 现象 ，Cavecchi 和 Spitkovsky ^ f£ 2020 年 模拟 了 爆发 过 程 中 
不 稳定 火焰 前 段 的 3D 模型 ， 不 稳定 的 火焰 前 端 导 致 整体 燃烧 速度 加 快 并 出 现 明显 火焰 涡 
流 ， 这 些 涡流 也 与 振荡 频率 的 漂移 有 关 ， 它 们 很 可 能 是 激发 表面 振荡 模式 的 原因 。 


模型 1， 火 焰 停 止 莹 延 


b) 
模型 2， 激 发 大 尺度 波浪 


» 


火焰 从 点 火 位 置 开 始 墓 延 


a) 


c) 
ik: 爆发 振荡 需要 存在 温度 不 对 称 性 ， 通 过 恒星 的 旋转 调制 ， 形 成 接近 自转 频率 的 脉冲 。 图 a) 表示 点 火 从 一 点 
始 ， 随 后 火焰 蔓延 至 整个 中 子 星 ， 其 中 红色 代表 高 温 ， 蓝 色 代表 低温 。 之 后 发 生 的 物理 过 程 目前 还 不 明确 ， 其 
中 一 种 可 能 性 是 图 b) 燃烧 物质 熄灭 ， 只 有 一 部 分 燃烧 物质 被 燃烧 ， 留 下 热 斑 一直 持 续 到 暴 尾部 ， 另 一 种 是 图 c) 
火焰 在 中 子 星 表面 扩散 ， 并 在 燃烧 物质 中 激发 起 大 尺度 的 波浪 。 


Els ”爆发 振荡 简化 示意 图 ” 
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现在 一 些 暴 的 振荡 信号 还 很 难 观测 到 ， 未 来 由 中 国 主导 的 eXTP 卫星 (the enhanced 
X-ray timing and polarimetry mission), A EH K 射线 时 变 观测 能 力 ， 将 为 研究 暴 振 荡 
机 制 提供 更 多 的 条 件 与 机 会 。 
5.1 ADHRA 

为 解释 暴 振荡 频率 向 上 漂移 和 振荡 主要 出 现在 暴 尾 部 的 现象 ， 人 们 提出 了 两 种 热 斑 模 
型 : 热 斑 传 播 即 热 斑 在 暴 上 升 阶段 中 临时 产生 ， 火 焰 从 点 火 点 扩散 到 整个 中 子 星 ， 另 一 种 是 
将 燃烧 限制 在 中 子 星 表面 的 一 小 片区 域内 ， 热 斑 持 续 存在 。 

Strohmayer 等 人 ”首先 提出 ， 火焰 在 中 子 星 上 扩散 时 的 水 平 运动 可 能 会 导致 振荡 频率 
的 漂移 ， 为 进一步 解释 频率 总 是 逐渐 上 升 的 现象 ， 提 出 另 一 种 假设 ; 在 加 热 燃 烧 层 的 垂直 脱 
胀 过 程 中 保持 角 动 量 守恒 ， 快 速 膨 胀 会 使 热 斑 旋 转 得 更 慢 ， 从 而 使 观测 到 的 频率 低 于 自转 速 
率 ， 并 且 随 着 层 的 冷却 和 收缩 频率 会 升 高 ， 逐 渐 达 到 自转 频率 。 在 大 多 数 暴 尾部 ， 拟 合 暴 能 
谱 所 获得 的 黑体 半径 与 中 子 星 的 半径 相似 ， 这 说 明火 焰 已 经 在 整个 中 子 星 表面 燃烧 ”。 但 
是 积累 爆发 物质 需要 数 小 时 甚至 数 天 ， 而 热 核 失控 不 到 1 s”， 为 了 同时 满足 整个 中 子 星 表 
面 的 点 火 条 件 ， 热 状态 必须 保持 一 致 ”， 而 吸 积 中 存在 一 些 不 对 称 性 现象 (例如 由 于 磁 通道 
或 赤道 边界 层 的 存在 )， 使 得 这 种 情况 (尤其 是 在 AXMP F) 不 太 可 能 发 生 。 所 以 ， 如 果 火 
焰 在 某 一 点 开始 燃烧 ， 那 么 火焰 需要 在 整个 中 子 星 表面 扩散 ， 但 是 目前 还 难以 解释 燃烧 物质 
在 这 时 可 以 达到 所 需要 的 热 均衡 状态 。 
散布 热 斑 模型 随 着 整个 暴 能 谱 黑 体 半径 增加 ， 热 斑 尺 度 增 大 ， 振 荡 幅 度 下 降 ， 
可 以 解释 暴 的 上 升 阶 段 中 存在 爆发 振荡 ， 但 是 无 法 解释 大 多 数 暴 振荡 信号 都 出 现在 暴 的 尾 
部 。 其 中 一 个 解释 是 因为 观测 效应 ， 暴 上 升 阶段 持续 时 间 很 短 ， 很 难 检测 ”。XBO 探测 难 
易 程 度 的 另 一 个 因素 是 点 火 处 的 纬度 ， 在 快速 旋转 的 中 子 星 上 ， 均 道 处 的 有 效 重力 较 低 ， 吸 
积 率 最 大 ， 这 样 更 容易 积累 点 火 所 需要 的 物质 ”。 研 究 发 现 ， 赤 道 处 的 I 型 X 射线 景 ， 会 
更 容易 探测 到 爆发 振荡 ”。 但 是 ， 在 边界 层 附近 ， 吸 积 物 质 的 核反应 过 程 介 于 稳定 燃烧 与 
不 稳定 燃烧 之 间 ， 可 能 在 高 纬度 处 点 火 ，”。 对 磁场 较 高 的 源 ， 由 于 局 部 密度 过 大 或 者 过 度 
m m S mn TN E 
焰 必 须 外 e de 衬 传 播 。 多 年 以 来 ， 控 制 火焰 扩散 的 过 程 在 I 型 X 射线 暴 中 一 
直 是 一 个 悬而未决 的 问题 ， 它 受热 传导 、 洲 流 和 对 流 等 各 种 传 热机 制 的 影响 "。Spitkovsky 
ey 特别 是 科 里 奥 利 力 与 热 燃烧 物质 膨胀 之 间 
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扩散 非常 重要 ， 因 为 热 物 质 的 快速 膨胀 和 横向 扩散 可 能 会 导致 火焰 失速 ”。 但 是 ， 科 里 奥 
利 力 约 束 对 于 ICR J17480 2446 (自转 频率 仅 为 11 Hz) EAE. 

同样 暴 尾部 的 高 振幅 振荡 现象 也 没有 准确 的 解释 ， 第 一 种 机 制 是 热 斑 得 以 保留 ， 燃 烧火 
焰 并 未 扩散 到 整个 中 子 星 上 。 让 热 斑 保 留 下 来 的 一 种 方法 是 通过 中 子 星 磁场 限制 燃料 。 吸 积 
物质 沿 磁力 线 进入 中 子 星 磁 极 ， 压 力 足 够 大 时 场 线 变形 ， 吸 积 物质 才能 越过 磁力 线 。Brown 
和 Bildsten"" 的 结果 表明 ， 要 使 物质 约束 在 正常 1 型 久 射 线 暴 的 点 火 深度 内 ， 其 场 强 至 少 应 
为 106 T， 这 要 远 远 高 于 LMXB 的 磁场 强度 ””， 而 磁 流 体 动力 学 的 不 稳定 性 也 可 能 使 磁 
约束 无 效 ” ， 因 此 ， 这 种 机 制 的 可 能 性 并 不 大 。 另 一 种 机 制 是 限制 火焰 前 端 ， 燃 烧 一 旦 开 
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台 ， 草 延 一 段 距离 后 停止 传播 ， 这 需要 以 某 种 方式 阻止 其 热 传 递 。Cavecchi 等 人 ™ 认 为 燃 
烧 物 质 与 磁场 的 径 向 分 量 相互 作用 ， 充 分 放大 磁场 的 水 平分 量 后 ， 也 能 阻止 热传导 。 对 于 低 
磁场 的 LMXB， 这 些 过 程 难以 出 现 ， 除 非 在 爆发 过 程 中 有 某 些 机 制 (例如 对 流 发 电 ) 暂时 增 
加 磁场 ” 。 
5.2 ”表面 模式 振荡 模型 

表面 模式 振荡 模型 认为 ,，I 型 X 射线 暴 点 火 以 及 火焰 在 中 子 星 周围 的 扩散 可 能 会 激发 中 
子 星海 洋 (由 离子 和 电子 的 等 离子 体 组 成 ， 分 布 在 中 子 星 外 层 到 等 离子 体 可 以 结晶 成 固体 的 
深度 ) 中 的 大 尺度 波浪 ， 海 洋 中 不 同 高 度 处 温度 有 差异 ， 以 致 X 射线 的 亮度 不 同 ， 通 过 星体 
转动 ， 表 现 为 XBO。 人 们 借鉴 了 地 球 海洋 表面 模式 振荡 的 研究 方法 ， 写 出 了 旋转 球体 振荡 
的 流体 动力 学 连续 性 方程 (质量 守恒 )， 动 量 和 能 量 守 恒 方 程 ， 并 简化 为 拉 普 拉 斯 潮汐 方程 
式 ”™"。 这 个 系统 的 恢复 力 是 浮力 和 科 里 奥 利 力 。 拉 普 拉 斯 潮汐 方程 的 非 轴 对 称 模式 解 依 
HF exp(imq) 的 方位 角 p， 其 中 方位 角 特 征 值 m 是 一 个 整数 ， 给 定 一 个 中 子 星 的 旋转 模式 
具有 方位 角 数 量 m 和 频率 v, 的 模式 ， 观 测 到 的 频率 为 : 


Vo = MVs FVr ， (14) 


其 中 ，z 是 中 子 星 自转 频率 ， 而 六 的 符号 是 正 号 还 是 负 号 ， 取 决 于 模式 是 沿 自转 方向 (向 
东 ) 还 是 道 自转 方向 (向 西 ) 传播 。 如 果 v. 取决 于 海洋 的 热 状 态 ， 就 可 以 解释 XBO 频率 漂移 
现象 ， 即 振荡 频率 会 随 着 爆发 的 进行 和 尾部 的 冷却 而 改变 。 

拉 普 拉 斯 潮汐 方程 式 有 三 种 模式 。 

(1) 庞 加 莱 模 式 。 向 逆 方 向 (向 西 ) 传播 ， 并 在 不 旋转 的 临界 处 变 为 纯 重 力 波动 ( 仅 靠 浮 
力 恢 复 )， 因 此 通常 被 称 为 g Bc 。 

(2) 开尔文 模式 。 沿 与 自转 相同 方向 (向 东 ) 传播 ， 在 不 旋转 的 临界 处 变 为 纯 重 力 波动 。 
但 与 庞 加 莱 模 式 不 同 的 是 ， 开 尔 文 模式 处 于 地 转 平衡 状态 (水 平 气压 梯度 力 与 科 里 奥 利 力 相 
平衡 ) ， 因 此 只 涉及 到 方位 角 运 动 ， 并 只 限于 在 赤道 上 进行 火焰 传播 。 

(3) Rossby 模式 。 沿 逆 方 向 (向 西 ) 传播 ， 也 称 为 + 模式 ”， 在 零 浮 力 的 极限 下 ，Rossby 
模式 是 由 科 里 奥 利 力 的 纬度 变化 驱动 的 。 在 不 旋转 的 临界 处 变 为 零 频率 的 一 般 解 。 如 有 不 同 
燃烧 物质 之 间 的 浮力 存在 ， 则 也 将 起 到 恢复 模式 的 作用 ， 因 此 也 通常 称 为 浮力 r 模式 。 

Heyl” 认 为 浮力 + 模式 是 最 有 可 能 的 ， 例 如 : 逆向 传播 的 模式 中 e| 仅 有 几 赫 兹 ， 给 
出 的 观测 频率 略 低 于 自转 频率 ， 随 着 燃烧 层 冷却 ，|v,| 变 小 ， 观 测 到 的 频率 将 趋向 于 自转 
频率 。 而 庞 加 莱 模 式 (或 g 模式 ) 预计 的 |v| 要 大 很 多 ， 所 以 被 排除 "”。 开 尔 文 模式 被 排 
除 原因 是 频率 将 随 着 海洋 冷却 而 下 降 ， 这 与 大 多 数 观测 到 的 频率 漂移 相反 ， 开 尔 文 模式 和 
庞 加 莱 模 式 也 都 将 燃烧 严格 限制 在 快速 旋转 的 中 子 星 上 的 赤道 区 域 ， 这 将 降低 整体 暴 振荡 
幅度 。 浮 力 r 模式 预言 燃烧 层 冷 却 时 的 频率 漂移 约 10 Hz， 比 观测 值 高 一 个 量 级 。Piro 和 
Bildsten 提出 振荡 出 现在 燃烧 层 底部 与 中 子 星 壳 层 之 间 的 区 域 ， 振 荡 频 率 漂 移 约 4 Hz， 依 
然 高 于 AMXP 源 的 测量 值 。Heyl upciteHey04 指 出 浮力 模式 如 果 限 制 在 光 球 层 ， 则 漂移 
频率 可 能 更 小 ， 细 节 有 竺 研究。 一 些 暴 的 光度 会 暂时 超过 爱 丁 顿 光度 ， 热 能 转换 为 动能 和 
引力 势能 ， 使 中 子 星 最 外 层 发 生 膨 胀 ， 称 为 光 球 半径 膨胀 (photospheric radius expansion, 
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PRE)， 通 常会 从 PRE 暴 中 检测 到 爆发 振荡 信号 ， 但 是 在 PRE 阶段 本 身 很 少 检测 到 爆发 振 
荡 。 为 了 解释 这 种 在 峰值 阶段 没有 检测 出 振荡 信号 的 现象 ，Cooper” 考虑 燃烧 层 高 效 对 流 
能 量 传输 机 制 ， 分 析 了 单 区 域 模型 ， 发 现 对流 确 实 抑制 了 振荡 的 发 生 ， 可 以 解释 一 些 暴 振荡 
TE AE Fe SE BO ILS 
Chambers 和 Watts ^ 考虑 浮力 r Bi] 2:5 zi dt S USURIS] AR, VETERA 
频率 漂移 约 1~3 Hz， 与 观测 结果 相符 ， 但 是 给 出 的 TBO 频率 与 自转 频率 的 偏差 仍然 比 观 
测 值 大 。 
5.3 ”冷却 尾 流 (cooling wake) 模型 
目前 在 一 些 X 射线 暴 的 尾部 检测 到 振幅 很 大 的 爆发 振荡 信号 ， 因为 扎 少 明显 的 高 梯度 
温度 差 ， 普 通 的 热 斑 模 型 无 法 解释 这 种 现象 。Cumming 和 Bildsten” 认 为 可 能 是 由 于 中 子 
星 表面 冷却 而 形成 温度 梯度 导致 了 暴 尾部 的 振荡 现象 ， 这 即 为 冷却 尾 流 模型 。 
Mahmoodifar 和 Strohmayer” 提出 了 三 种 冷却 尾 流 模型 ， 其 中 “规范 ”冷却 尾 流 模型 
分 为 两 种 : 一 种 假设 中 子 星 上 每 一 片区 域 以 相同 的 方式 燃烧 和 冷却 ， 中 子 星 表面 一 旦 从 某 点 
开始 点 燃 燃烧 ， 火 焰 扩 散 到 每 个 区 域 之 后 该 区 域 温度 会 由 指数 上 升 和 下 降 ， 使 得 中 子 星 表面 
出 现 一 定 的 温度 梯度 ， 温 度 梯 度 和 中 子 星 自 旋 相 影响 ， 使 得 暴 尾部 出 现 爆 发 振荡 信号 。 第 二 
种 “规范 ”模型 称 为 纬度 依赖 模型 ， 因 为 该 模型 假设 中 子 星 表面 的 冷却 时 间 随 纬度 变化 而 变 
化 ， 在 快速 旋转 的 中 子 星 表 面 ， 有 效 重力 加 速度 沿 赤 道 方 向 减 小 ， 在 极点 处 最 大 ， 所 以 质量 
吸 积 率 和 冷却 时 间 在 赤道 处 最 大 、 极 点 处 最 小 。 虽 然 这 个 模型 考虑 了 冷却 时 间 变 化 因素 ， 但 
保持 了 自 旋 轴 对 称 性 。 第 三 种 为 非 对 称 冷却 唯 象 模型 ， 冷 却 时 间 变 化 不 再 满足 中 子 星 自 旋 
对 称 性 ， 也 不 考虑 每 个 单独 区 域 的 温度 变化 ， 而 是 只 考虑 已 燃烧 区 域 ( 热 区 ) 和 剩 下 区 域 
( 冷 区 ) 的 温度 。 火 焰 燃 烧 到 整个 中 子 星 表面 时 ， 时 间 很 短 (用 时 大 概 为 1 s)， 中 子 星 表 面 温 
度 梯 度 可 忽略 ， 当 开始 冷却 时 ， 先 燃烧 的 区 域 最 先 冷却 ， 假 设 冷却 与 燃烧 以 相同 方式 扩散 ， 
但 是 时 间 更 长 (用 时 大 概 为 10 s)， 因 此 先 燃 烧 区 域 不 仅 先 冷却 ， 而 且 还 要 比 其 他 区 域 冷却 时 
间 短 很 多 ， 此 过 程 可 以 增加 暴 尾部 振荡 信号 的 振幅 。 
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6 其 振荡 的 应 用 


€ 6.1 ”测量 中 子 星 的 质量 和 半径 


精确 测量 中 子 星 的 质量 和 半径 对 于 我 们 认识 致密 物质 状态 方程 和 中 子 星 的 内 部 结构 非 
常 关键 。LMXB 的 中 子 星 质量 和 半径 可 由 多 种 独立 的 方法 测定 ， 例 如 ， 由 光 球 层 膨胀 暴 的 
时 间 分 辩 能 谱 号 ， 测 量 质量 和 半径 ， 已 知 距离 ， 拟 合 宁静 态 的 热 谱 瑟 ， 测 量 质量 和 半径 ; 
中 子 星 表 面 特征 谱 线 ”测量 质量 和 半径 之 比 等 。 拟 合 脉 冲 轮 廓 也 能 测量 中 子 星 的 质量 和 
半径 。X 射线 脉冲 星 因 表面 热 斑 的 周期 性 自转 ， 表 现 出 脉冲 形式 的 辐射 。 在 强 引 力 场 的 影响 
下 ， 热 斑 的 辐射 束 弯曲 ， 因 此 观测 到 的 脉冲 轮 廊 ， 与 脉冲 星 的 质量 与 半径 之 比 有 关 ， 同 时 也 
与 热 斑 形状 、 与 视线 方向 的 夹 角 等 因素 相关 。Psaltis 等 人 四 研究 了 快速 转动 的 吸 积 X 射线 
脉冲 星 (自转 频率 180 ~ 620 Hz)， 此 时 星体 因 高 速 旋 转 而 形变 ， 明 显 改 变 了 脉冲 轮廓 ， 这 
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个 效应 独立 于 光线 弯曲 ， 给 出 了 另外 的 质量 半径 的 限制 关系 ， 从 而 可 以 不 用 测量 源 的 距离 
而 测量 出 星体 的 质量 和 半径 。NICER 卫星 的 主要 科学 目标 就 是 研究 转动 供 能 X 射线 脉冲 星 
(脉冲 星 的 辐射 能 量 来 源 于 星体 的 自转 能 ) 的 物 态 ， 目 前 ，NICER 累积 了 大 量 的 观测 数据 ， 
测量 了 PSR J0030+0451 的 质量 和 半径 ”， 约 束 了 致密 物质 的 状态 ， 此 方法 依赖 于 辐射 光 
子 来 自 于 中 子 星 表面 而 不 是 中 子 星 附近 的 磁 层 。 在 X 射线 脉冲 星 中 应 用 的 脉冲 轮廓 拟 合 方 
法 也 适用 于 XBO， 高 精度 质量 和 半径 的 测量 ， 需 要 探测 到 106 个 光子 。 未 来 eXTP 卫星 将 
会 搭载 大 有 效 面积 的 探测 器 ， 通 过 测量 中 子 星 LMXB 的 XBO， 累 积 多 次 观测 得 到 的 脉冲 轮 
廊 ， 拟 合 其 脉冲 轮廓 ， 将 对 中 子 星 物 态 给 出 更 加 严格 的 限制 。 现 在 已 经 观测 到 了 30 颗 中 子 
Æ LMXB ( 吸 积 供 能 ) 的 XBO 现象 ， 观 测 样本 数 远 远 多 于 转动 供 能 X 射 线 脉冲 星 ， 更 多 的 
中 子 星 质量 和 半径 的 测量 将 在 更 宽 的 质量 谱 上 约束 致密 物质 状态 名。 结合 其 他 方法 测量 到 的 
星体 质量 和 半径 ， 将 有 效 地 减 小 质量 和 半径 测量 的 误差 。 
6.2 ”测量 中 子 星 自转 

多 年 来 ， 人 们 通过 多 波段 观测 搜索 各 类 中 子 星 ， 研 究 中 子 星 的 自转 频率 分 布 和 演化 ， 因 
为 它们 是 理解 恒星 和 双星 系统 演化 的 重要 因素 。 观 测 表明 ， 中 子 星 LMXB 的 XBO 频率 接 
近 于 星体 的 自转 频率 ， 其 误差 在 几 赫 效 之 内 ， 由 XBO 的 观测 ， 独 立 发 现 了 较 多 吸 积 供 能 
子 星 的 自转 频率 。 通 过 对 中 子 星 LMXB 4U 1608-522 的 XBO 观测 ， 发 现 其 自转 频率 约 为 
620 Hz， 这 是 由 X 射线 能 段 测量 到 转动 最 快 的 中 子 星 。 发 现 更 多 快速 自转 的 中 子 星 ， 对 理 
解 射 电 毫 秒 脉冲 星 的 再 循环 机 制 至 关 重 要 ， 尤 其 是 吸 积 演化 过 程 中 的 物质 转移 、 磁 场 衰减 和 
吸 积 加速 等 过 程 。 另 外 一 个 重要 的 问题 是 中 子 星 可 以 达到 的 最 大 自转 频率 四， 亚 毫秒 周期 的 
中 子 星 将 严格 而 干净 地 约束 致密 物质 的 状态 方程 ”。 而 且 ， 在 快速 旋转 中 子 星 上 ， 是 否 
会 观测 到 暴 振荡 ， 这 也 是 目前 尚未 解决 的 问题 。 


7 结论 与 展望 


本 文 介绍 了 中 子 星 LMXB 中 的 XBO 现象 ”人 们 提出 了 多 个 理论 模型 解释 XBO 的 观测 
现象 ， 但 这 些 模 型 难以 解释 所 有 的 观测 现象 ， 尚 待 完善。 目前 ，NICER 探测 到 了 XBO I 
象 ， 并 将 观测 到 更 多 的 I 型 X 射线 暴 。 未 来 ， 中 国 主导 的 eXTP 卫星 也 将 观测 中 子 星 作为 
主要 的 研究 目标 ， 我 们 期 待 通过 观测 XBO， 发 现 亚 毫 秒 脉冲 星 ， 精 确 测量 中 子 星 的 质量 和 
半径 并 严格 限制 致密 物质 状态 ， 建 立 更 加 完善 的 XBO 理论 模型 。 


?例如 ， 可 能 发 现 大 质量 (> 2M) 或 小 质量 (< 1Mo) 的 中 子 星 ， 而 小 质量 中 子 星 的 半径 如 能 准确 测量 ， 将 有 效 地 区 分 奇 
F 星 和 中 子 星 模型 。 


9 是 否 存 在 亚 毫秒 中 子 星 
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Research Progress of X-ray Burst Oscillation 


LI Qing-xin, LI Zhao-sheng 
(Key Laboratory of Stars and Interstellar Medium, Xiangtan University, Xiangtan 411105, China) 


Abstract: X-ray burst oscillations are coherent signals which are detected in neutron star 
low-mass X-ray binaries during their type-I X-ray bursts. The frequency of X-ray burst 


oscillation is close to the neutron star spin rate. In recent years, about 30 neutron star 
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low-mass X-ray binaries have been detected X-ray burst oscillations from RXTE, Swift and 
NICER observations. They are important to probe the dense matter equation of state and 
the effects of strong gravity, to study the high magnetic field environment on the neutron 
star surface. The studies of X-ray burst oscillations are also crucial to understand the spin- 
up process of neutron star low-mass X-ray binaries through accretion and the evolution into 
millisecond pulsars. The different characteristics of burst oscillation signals between persis- 
tent accretion-powered pulsars, intermittent accretion-powered pulsars and sources without 
detectable accretion-powered pulsations can help us to study more about it. We briefly in- 
troduce type I X-ray bursts, and then specifically explain the methods to search for burst 
oscillations and provide the observational results of X-ray burst oscillations. The theoretical 
burst oscillation model, i.e., the hot spot model, the surface mode model and the cooling 
wake model, are discussed. At present, we note that these models are difficult to explain 
all the observed phenomena of X-ray burst oscillations. In future, the China-led mission, 
the enhanced X-ray Timing and Polarimetry (eXTP) observatory will observe neutron star 
low-mass X-ray binaries and search X-ray burst oscillations. By modeling the stacking pulse 
profile of X-ray burst oscillation, the mass and radius of neutron star will be measured ac- 
curately, and the neutron star equation of state will be tightly constrained. We also expect 
more neutron stars will be measured their spin frequency, if luck, sub-millisecond pulsar will 


be discovered. 


Key words: X-ray burst; X-ray burst oscillation; Low-mass X-ray binary; Pulsar 


附录 A 


搜索 Swift/BAT B X S] 2c 25d E] X 射线 暴 振荡 信号 ， 以 Swift 在 2006 年 观测 到 的 
4U 0614+09 的 爆发 数据 为 例 ， 步 又 如 下 。 

(1) 下 载 对 应 观测 数据 。 

(2) 在 数据 所 在 目录 ， 终 端 输入 命令 


batgrbproduct indir=00234849000 outdir=../00234849000-results 


对 下 载 数据 进行 处 理 ， 输 出 文件 按 类 型 排放 在 输出 目录 中 。 
(3) 在 目录 00234849000-results/events， 得 到 能 量 范围 为 13 ~ 20 keV 的 事例 数据 转变 
为 时 间 分 辩 率 为 0.25 s 的 并 道 数据 ， 终 端 输入 : 


batbinevt infile-sw00234849000b all.evt outfile-rate no. 13. 20keV O25s.1c 


outtype-LC  timedel-0.25 timebinalg-u  energybins-13-20 detmask= 
../auxil/sw00234849000b qmap.fits ^ ecol-ENERGY . weighted-YES 


208 天 文学 进展 39 3$ 


outunits-" RATE?! clobber-yes 


良 像素 图 ) 文件 名 称 ， 用 来 反映 BAT 望远镜 中 科学 数据 


其 中 detmask 是 探测 器 的 质量 图 
质量 的 任何 操作 问题 。 

(4) 画 出 上 面 得 到 并 道 数据 的 光 变 曲线 ， 查 看 X 射线 暴 所 在 时 间 范 围 ， 爆 发 时 间 内 光 
变 曲线 取 每 4 s 进行 快速 傅 里 叶 变 换 ， 每 次 移动 步 长 为 0.25 s， 得 到 功率 密度 谱 ， 本 例 用 
Python 对 数据 进行 相关 操作 ， 有 具体 代码 见 代 码 块 (1)。 


E 


Listing 1  FFT.py 


from astropy.io import fits 

import numpy as np 

import matplotlib.pyplot as plt 
from stingray import Powerspectrum 


from stingray.events import EventList 


data = fits. open( file) #tile 为 事例 数据 文件 

ev time = np.array(data[1].data.field("TIME")) # 光 子 到 达 时 间 
energy = data[1].data.field("ENERGY") # 光 子 对 应 能 量 范 转 
en index = np.where((energy >= 13) & (energy «- 20)) 

ev en time = ev time[en index] # 取 出 13-20KeV 能 量 范围 的 光子 


time = ev en time - ev_en_time[0] 
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ev = EventList (time) 

lc = ev.to lc(dt-0.1) 

t_start = 0 

t end = time[-1] 

t shift = 0.25 # 时 间 移 动 步 长 


length = int((t end - t start) / t shift) + 1 


[o E 
O 0 - 0 O0 w^ 


p = np.zeros(length) 


N 
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t mark - t start 


ND 
pà 


t_medium = np.zeros(length) 
t_bin = 4 # 时 间 间 隔 


N N N 
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for i in range(length): 
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index = np.where((time >= t_mark) & (time <= t_mark + t_bin)) 


N 
[o>] 


t medium[i] = t mark + t bin / 2 


N 
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times - time[index] 
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ev = EventList(times) 
lc new = ev.to lc(dt-1e-4) 
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ps_leahy = Powerspectrum(lc new, norm-'leahy') 


w 
nes 


plt.rc('font', family-'Times New Roman', size-12) 
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.plot(ps leahy.freq, ps leahy.power, lw-3) 
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.tick params(axis-"x", direction-"in", labelsize-10) 
.tick params(axis-"y", direction-"in", labelsize-10) 
.xlim(11, 900, 0.2) 

plt.ylim(0, 100) 


t mark = t mark + t shift 
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(5) 得 到 功率 密度 谱 在 414 Hz 附近 有 明显 的 振荡 信号 ， 取 频率 范围 413 ~ 416 Hz，13 
~ 20 keV 的 事例 数据 ， 做 22 检验 ， 将 Z2 —22, 26, 30, 34, 38 的 等 高 线 投 影 在 光 变 曲线 上 ， 


得 到 的 动态 功率 谱 与 Strohmayer 等 人 四 对 4U 0614+09 爆发 数据 的 处 型 


P RES MEE! 


号 频率 为 414.75 Hz， 具 体操 作 见 代码 块 (2)。 


I f 


Listing 2  Zn.py 
1 |from astropy.io import fits 
2 |import numpy as np 
3 | import matplotlib.pyplot as plt 
4 |from stingray import Lightcurve 
5 |from stingray.pulse.search import z n search 
6 |from stingray.events import EventList 
7 |from matplotlib.pyplot import MultipleLocator 
8 
9 |data = fits. open(ev file) # 打 开 事 例 数据 文件 
10 |ev time = data[1].data.field("TIME") 
11 |energy = data[1].data.field("ENERGY") 
12 |en index = np.where((energy >= 13) & (energy <= 20)) 
13 |ev en time - ev time[en index] 
14 |time = ev en time - ev en time[0] 
15 |t start = 0 
16 |t end = time[-1] 
17 |t shift = 0.25 
18 |length = int((t end - t start) / t shift) + 1 
19 |frequencies = np.arange(413, 416, 0.1) # 搜 索 频 率 范 围 
20 |zsquare = np.zeros(( len(frequencies), length)) 
21 |t mark = t start 
22 |t medium - np.zeros(length) 
23 |x zn = MultipleLocator(1) 
24 |y zn = Multiplelocator(1) 
25 |t bin- 4 
26 
27 |for i in range(length): 
28 index = np.where((time >= t mark) & (time <= t mark + t bin)) 
29 t medium[i] = 七 mark + t bin / 2 
30 times - time[index] 
31 ev = Eventlist(times) 
32 t mark = t mark + t shift 
33 nharm - 1 
34 freq, zstat = z n search(ev.time, frequencies, nharm = nharm) 
35 plt.plot(freq, zstat) 
36 zsquare[:, i] = zstat 
37 
38 |lc 1 data = fits. open(file lc 025) # 打 开 并 道 数据 文件 


202306.00346v1 


chinaXiv 


210 


39 
40 
41 
42 
43 
44 
45 
46 
4T 
48 
49 
50 
51 
52 
53 
54 
55 
56 
57 
58 


time lc 1 = 1c_1_data[1] .data.field("TIME") 

counts lc 1 = lc 1 data[1].data.field("RATE") 

time lc 1 = time lc 1 - time lc 1[0] 

lc p = Lightcurve(time lc 1, counts lc 1) 

lc p.plotO 

plt.rc('font', family-'Times New Roman', size-12) 

fig, axi = plt.subplots() 

axi.tick params(axis-"x", direction-"in", labelsize-10) 
axi.tick params(axis-"y", direction-"in", labelsize-10) 
axi.plot(lc p.time, lc p.counts, 'k', lw-1) 

axi.set xlabel('Time (s)') 

axi.set ylabel('Counts/s/cm$^2$') 

ax2 = axi.twinx() 

levels = [22, 26, 30, 34, 38] # 设 置 等 高 线 参 数 

X, Y = np.meshgrid(t medium, frequencies) 

ax2.tick params(axis-"y", direction-"in", labelsize-10) 
ax2.contour(X, Y, zsquare, levels) 

ax2.set ylabel('Frequency(Hz)') 

ax2.set ylim(413, 416) 

fig.tight layoutO 
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